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ABSTRACT 

The optical-UV continuum of quasars is broadly consistent with the emission 
from a geometrically thin optically thick accretion disk (AD). The AD produces the 
ionizing continuum which powers the broad and narrow emission lines. The maximum 
AD effective temperature is given by Toff, max = fmnxiM /Ad'^Y^'^ , where M is the black 
hole mass, M the accretion rate, and /max is set by the black hole spin a^,. For a low 
enough value of M /M'^ the AD may become too cold to produce ionizing photons. 
Such an object will form a lineless quasar. This occurs for a local blackbody (BB) 
AD with a luminosity Lopt = lO"**^ erg s"^ for M > 3.6 x WMq, when a* = 0, 
and for M > 1.4 x 10^°Mq, when a* = 0.998. Using the AD based M, derived 
from M and Lopt, and the reverberation based M, derived from Lopt and the H/3 
FWHM, V, gives Teff^max oc L~p°'^^D~^-^^. Thus, Tcff^max is mostly set by v. Quasars 
with a local BB AD become lineless for v > 8,000 km s^^, when a* = 0, and for 
V > 16,000 km s"""^, when a* = 0.998. Higher values of v are required if the AD is 
hotter than a local BB. The AD becoming non- ionizing may explain why line emitting 
quasars with v > 10, 000 km s~^ are rare. Weak low ionization lines may still be present 
if the X-ray continuum is luminous enough, and such objects may form a population 
of weak emission line quasars (WLQ). If correct, such WLQ should show a steeply 
falhng SED at A < lOOOA. Such an SED was observed by Hryniewicz et al. in SDSS 
J094533.99+100950.1, a WLQ observed down to 570A, which is well modeled by a 
rather cold AD SED. UV spectroscopy of z 1 — 2 quasars is required to eliminate 
potential intervening Lyman limit absorption by the intergalactic medium (IGM), and 
to explore if the SEDs of lineless quasars and some additional WLQ are also well fit 
by a cold AD SED. 

Key words: accretion, accretion disks — black hole physics — galaxies: active — 
galaxies: quasars: general 



1 INTRODUCTION 

The optical-UV spectral energy distribution (SED) of 
quasars is broadly consistent with the expected emission 
of a thin accretion disk (AD) around a massive black hole 
(Shields 1978; Czerny & Elvis 1987; Sun & Malkan 1989; 
Laor 1990; Blaes et al. 2001; Shang et al. 2005). Various 
spectral features, such as the general turnover at A < lOOOA 
(Zheng et al. 1997; Teller 2002; Barger & Cowie 2010), the 
observed optical-UV spectral slopes (Bonning et al. 2007; 
Davis et al. 2007), the small dispersion in the UV/optical 
flux ratio (Davis & Laor 2011), and microlensing variabil- 
ity (Morgan et al. 2010; Blackburne et al. 2011), may imply 
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various modifications beyond the simplest AD models (see 
review by Koratkar & Blaes 1999, and citations thereafter). 
One should note that some of the expected AD spectral 
features may be diluted by various unrelated emission com- 
ponents (Kishimoto et al. 2004; 2008). 

The peak of the simplest AD model of local blackbody 
(BB) emission occurs at fc'peak (X (M/M^)^/* (see §2), or 
equivalently z/peak oc {l/M)^^*, where M is the black hole 
mass, Al is the accretion rate, and / = L/I/Edd is the ratio 
of the bolometric luminosity to the Eddington luminosity. 
Thus, when M increases from WMq to 10^ M©, at a fixed 
I, the disk cools by a factor of 100. This factor is consistent 
with the drop from J/peak — 1—2 keV in X-ray binary systems 
(e.g. RemiUard & McClintock 2006), to Up^^y, ~ 10 - 20 eV, 
observed in quasars. The quasar SED peaks at about the 
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H ionization threshold, and thus provides ample ionizing 
luminosity which powers the broad and narrow line emission. 

A further extrapolation of the expected AD SED to 
M > 10® M© and / < 0.1 leads to i^poak < 10 eV and thus to 
a non-ionizing AD SED (e.g. see Laor & Netzer 1989 fig.5, for 
a model where t^peak ^1 — 2 eV). Such models appeared in- 
consistent with the observations, as the SED of such quasars 
cannot power strong line emission. Such quasars will appear 
mostly as luminous continuum sources. Some line emission 
can still be produced through X-ray photoionization. How- 
ever, in luminous quasars the fraction of the bolometric lu- 
minosity in the x-ray drops well below 10% (Just et al. 2007), 
and the fraction of the bolometric power available for line 
emission thus drops from > 50% to < 10%., i.e. by a fac- 
tor of > 5. The total line emission should drop by a similar 
factor. 

Observational evidence that weak line quasars (WLQ) 
do exist started to accumulate since the discover}' of McDow- 
ell et al. (1995) that in PG 1407+265 all lines, except Ha, 
are exceptionally weak. A BL Lac origin could be clearly 
excluded based on its SED which overlaps well the mean 
SED of optically selected quasars (Elvis et al. 1994). Addi- 
tional WLQ were discovered in various other studies (Fan 
et al. 1999, 2006; Anderson et al. 2001; Hall et al. 2002, 
2004; Reimers et al. 2005; Ganguly et al. 2007; Leighly et 
al. 2007a; Hryniewicz et al. 2010), and were followed up by 
dedicated studies of larger samples of WLQ at a range of 
wavelengths (Shemmer ct al. 2009, 2010; Diamond-Stanic 
et al. 2009; Plotkin et al. 2010a, 2010b; Wu et al. 2011). 
The observed radio to X-ray SEDs of most WLQ is consis- 
tent with the mean SED of radio quiet quasars, and clearly 
excludes the BL Lac interpretation for the absence of line 
emission. Proposed explanations for WLQ invoked unusual 
Broad Line Legion (BLR) properties, such as a low covering 
factor, an anisotropic ionizing source, or BLR shielded from 
the ionizing source. 

A possibly related issue is that of 'true type 2' AGN, i.e. 
AGN which lack broad lines, but appear to be unobscurcd 
(Tran 2003; Shi et al. 2010). These objects differ from WLQ 
in two ways. First, they are defined by the lack of broad 
lines, and not by the lack of lines in general. However, once 
« > 1, the significant narrow lines are redshifted from optical 
spectroscopy, and one cannot separate a lineless quasar from 
a type 2 AGN. Second, the 'true type 2' AGN tend to reside 
at I < 10"^ (Shi et al. 2010; Trump et al. 2011), in contrast 
with WLQ which typically have i ~ 0.1 — 1 (e.g. Shemmer 
et al. 2010). 

Various suggestions were made to explain the lack of 
broad lines in unobscurcd type 2 AGN. Nicastro (2000), and 
Nicastro et al. (2003), suggested that the BLR is formed 
by accretion disk instabilities occurring in a critical radius 
where the disk changes from gas pressure dominated to ra- 
diation pressure dominated. This critical radius can become 
smaller than the innermost stable orbit, and then the BLR 
cannot form. Czerny et al. (2004) and Elitzur & Ho (2009) 
suggested that the transition of the AD from a cold geo- 
metrically thin flow, to a hot advection dominated accre- 
tion flow, eliminates the AD wind which launches the BLR. 
Elitzur & Shlosman (2006) suggested that the BLR is formed 
by a torus outflow, and below a certain accretion rate the 
outflow transforms into radio jets, and the BLR disappears. 
The various physical mechanisms suggested above to elimi- 



nate the BLR generally occur below some values of I, which 
may be a function of M. These ranges correspond to a cer- 
tain maximum line width, which may be function of L, be- 
yond which all AGN become 'true type 2' AGN. Laor (2003) 
pointed out that the BLR appears not to be detectable at 
FWHM > 25, 000 km s"\ which may either be a detection 
limit, or a maximal velocity where the BLR can exist. How- 
ever, no physical mechanism was put forward to explain why 
the BLR FWHM may be the primary parameter for the ex- 
istence of the BLR. As shown below, a non-ionizing AD may 
provide a simple explanation for the observed luminosity in- 
dependent upper limit on the BLR FWHM. 

In the following section we derive the accretion disk 
model parameters which produce a cold AD, i.e. an AD with 
a non-ionizing SED. In §3 we show that the SED of a WLQ, 
observed as fax as 570A, is fit surprisingly well by a cold 
AD model. In §4 we discuss additional implications, and the 
relevance to other observed properties of AGN. The main 
conclusions are summarized in S5. 



2 NON-IONIZING AD MODEL PARAMETERS 

Below we first derive the relation between Toff, max and the 
accretion disk model parameters, based on the standard 
Shakura & Sunyaev (1973, SS73) thin disk model. The flux 
emitted per unit area is 

F^|^^mM,a.) (1) 

where R is the radius, and f{R, M, o*) is a dimensionless fac- 
tor set by the inner boundary condition, and the relativistic 

effects (Novikov & Thorne 1973; Riffert & Hcrold 1995). It 
is convenient to use the dimensionless radius, r = R/Rg, 
where Rg = GM/c^, which gives 

^=8&W^(^'-)' (2) 

where /(r, a,) < 1. 

The local effective temperature is Teff = (F/o-)!/*, and 
it scales as 



Teff = To/(r,a,)^/V-^/^ 
where 

^ _ / 3c« \ Ml/4 



yS-KG'^a J M1/2-' 
or in more convenient units 
To = 8.6 X 10^Mi^''*m8 
where 

M = Ml Mq yr~^, 
and 

M = 10* mg Mq, 

or in cgs units 

M = 6.3x10^^ Ml grs"\ 

and 

M = 1.989 X lO^Vs gr. 



K 



(3) 



(4) 



(5) 
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We denote by /max the maximum value of f{r,a,), which 
sets the maximum disk temperature Tefi.max, assuming a 
local BB emission. The ratio rcfi, max/To is then a function 
of a, only. Ratios for representative values of a, are provided 
in Table 1. As a, increases, the innermost marginally stable 
disk radius Vn-is decreases, and Tcs^max/To increases. 

Note that /max is only a function of a* as we have as- 
sumed the standard AD structure, where the viscous torque 
increases inwards, reaches a maximum, and then drops until 
it vanishes at the innermost stable circular orbit (ISCO). 
In reality, there is likely to be some magnetic stress at or 
near the inner boundary and continuing emission inside the 
ISCO. Quantifying the level of additional emission from near 
or inside the ISCO is an active research topic (e.g Kulkarni 
et al. 2011, Noble et al. 2011). However, any additional emis- 
sion is not likely to cause significantly greater variation in 
/max than varying a, alone. For example. Noble et al. (2011) 
generally argue for a larger effect due to emission inside the 
ISCO than Kulkarni et al. (2011). But they find that an 
a, — would appear as an a, — 0.2 — 0.3 model if the emis- 
sion at the disk surface is locally everywhere a BB. (Note, 
however, that the radiation inside the ISCO may be far from 
thermodynamic equilibrium due to the drop in disk surface 
density, so local BB emission may not be a good approxima- 
tion there.) Hence, in this work we will focus on the simple 
AD models where /max is determined solely by a*. 

What is the frequency i^max where the accretion disk 
vLi, peaks? For a single temperature BB, it is simple to show 
that the peak occurs at ft/^max/fcTmax = 3.92. A local BB AD 
is a superposition of blackbodies with T ^ Tmax, weighted by 
the emitting surface area, and convolved by the relativistic 
effects. This superposition results in an integrated spectrum 
with a broader spectral energy distribution (SED) compared 
to a single temperature BB, which peaks at /li-'max/fcrmax < 
3.92. Table 1 lists the values of /i/^max/fcTmax, corresponding 
to an AD observed at /i = cos^ = 0.8, where 6 is the incli- 
nation angle of the AD to the line of sight (/x = 1 is face 
on). 

How low should I'max be for an AD to be called non- 
ionizing? Let us define a non- ionizing AD, or a cold AD, 
as an AD where the ionizing luminosity is 0.01 of the 
bolometric luminosity. We now look for fo.oi, the frequency 
above which the integrated luminosity is 1% of the bolomet- 
ric luminosity, i.e. 



L^dv = 0.01. 



If i/Q.oi < 3.29 X 10^^ Hz, then the AD is cold. The values of 
J^o.oiJ^max depends on the AD SED, which is set by a,. Table 
1 provides vo.oi/i^mmx values for different a* values, which all 
cluster at ~ 3.5. Thus, a non-ionizing AD peaks at z^max < 
9.4 X 10^* Hz, or Amax > 3200A. Since the typical x-ray 
luminosity observed is Lx < O.lLboi, the ionizing luminosity 
of a cold AD is ~ 10 smaller than L^, and the ionization 
will be essentially all from Lx- 

We can also define a weakly-ionizing AD as an AD 
where the ionizing luminosity is ^ O.lLboi. For such an ac- 
cretion disk the AD ionizing luminosity is ~ L^,. Such ob- 
jects will form a transition from UV dominated ionization 
to X-ray dominated ionization. Table 1 lists the values of 
J^o.i/i^max, which clustcr at ~ 2. Setting uo.i < 3.29 x 10^^ Hz 



for a weakly-ionizing AD implies that in such an AD z^max < 

1.65 X 10^^ Hz, or Amax > I8OOA. 

Using the calculated values of I'o.oi/z^max, 
/iz^max/fcTcff.max, aud Tcfj^max/To, for a giveu a,, as 
provided in Table 1, we derive the following upper limits 
for a cold AD: 

To s; 2.08 X 10^ K 

for a* = 0, and 

To s; 7.98 X 10* K 

for a, = 0.998. Inserting these limits into eq. 5, yields the 
conditions 



Ml < 3.42 X 10"^mg 
when a, — 0, and 



Ml < 7.41 X 10"^mg. 



(6) 



(7) 



when a* = 0.998. If mg is known, then Mi can be estimated 
from the AD model fit to the optical (5100A) luminosity, 
Lopt. Using eq. 8 in Davis & Laor (201l|3 

Ml = 3.5mg-''-«»Lip=t,45, (8) 
where Z/opt,45 = Lopt/lO*^, we obtain the lower limits on 



^opt. 



below which the AD is cold. 



Lopt,45 < 9.85 X 10-^m^-^^ 

for a* = 0, and 

Lopt,45 < 7.66 X W-'^ml^^ 



(9) 



(10) 



for a* = 0.998. Thus, even luminous quasars can become 
lineless, if mg is high enough. For example, for Lopt ~ 
10*® erg s~^, the AD becomes too cold to ionize when 
M > 3.6 X 10^ M© for a. = 0, or for M > 1.4 x W^°Mq 
for a* — 0.998. The corresponding value for I are < 0.22 and 



< 0.06, using the approximation Lb 



lOL, 



opt , 



i.e. Lbni = 



10*' erg s~*, and the implied relation I = Lopt,45/1.25mg. 
Low luminosity AGN need to have lower / values to have a 
cold AD. For example, AGN at Lopt ~ 10*^ erg s~^ need to 
have ; < 2.6 X 10"^ when a, = 0, and Z < 7 x 10"* when 
a* = 0.998, to become cold. 

The value of M is commonly estimated using the broad 
emission lines. Kaspi et al. (2000) gives 

mg = 1.5I/op'l^45 1^3000, (11) 

where U3ooo=H/3 FWHM/3000 km s"\ Inserting the Kaspi 
et al. relation to eq. 8 above, gives (Davis & Laor 2011, eq. 
9) 

Ml = 2.5L»p«t^45Woo', (12) 
Inserting these two expressions to eqs. 6-7 yields 
«30oo > 2.72L;7p°;45**, (13) 
when a* = 0, and 



^ Note that there is an error in the numerical factors in eqs. 5 and 
7 of Davis & Laor (2011). The numerical factor in eq. 5 should be 
40/7r2(6/iG2/5c2)i/3 rather than 160/it'-' (Giv^hG^c'^ /5)^^'^ . The 
numerical factor in eq. 7 should be 1.4 M© yr~^ rather than 
2.6 Mq yr"'^. However, the fitting function for M in their eq. 8, 
which corresponds to our eq. 8, is unaffected. 
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i;3ooo ^ 5.28L-p1f4f (14) 

when a« = 0.998. Thus, the criterion for a non-ionizing AD 
is set almost purely by the value of the H/3 FWHM, and 
is very weakly dependent on iopt- The AD becomes cold 
for H/3 FWHM> 8160 km s-\ when a, = 0, or for H/3 
FWHM> 15, 840 km s"\ when a* = 0.998. The quasar may 
become lineless when the H/3 FWHM is above these limits. 
These maximal values for the H/5 FWHM for an ionizing AD 
may correspond to the maximum values of the observable H/? 
FWHM, as quasars with broader lines may drop out from 
quasar surveys which are based on emission lines selection, 
due to their significantly weaker line emission. 

Figure 1 shows the range of SEDs for various AD 
models. All models have the same absolute accretion rate 
Ml = 1, and fj, = 0.8. The upper two panels show local BB 
accretion disk models. In both panels the SED gets softer 
as M increases, since !^max « M^^^^ (sec eq. 4). At a given 
M, the SED is softer for the a, — model, and i^max is 
a factor of 2.4 lower compared to the a, — 0.998 model 
(as can be deduced from Table 1). Note that although the 
a* = 0.998 AD is hotter than the a* = case by a factor of 
Teff,max(a, = 0.998)/Teff ,,nax(a, = 0) = 4.5, the SED of the 
a* = model has fti^max/fcTeff^max = 2.18 compared to only 
1.18 for the a, = 0.998 model, which makes the a« = 0.998 
model fmax only a factor of 2.4 higher than for the a* = 
model. Each curve is labeled with the ionization fraction 
/ion, which is defined as the fraction of the bolomctric lumi- 
nosity that is radiated at all wavelengths shortward of the 
the Lyman edge. 

There arc two effects which can make the observed spec- 
trum harder than shown in Fig.l for the local BB mod- 
els. These effects will increase the required minimal H/3 
FWHM for an AD to become cold. First, the AD SED is 
non isotropic. In the Newtonian case only the normaliza- 
tion of the AD luminosity changes with inclination. How- 
ever, the relativistic effects change also the shape of the 
SED, making it harder closer to an edge on view. The BLR 
must subtend an appreciable fraction of the sky, as seen by 
the ionizing source. The BLR is most likely in the form of 
a thick torus- like configuration, coplanar with the AD. The 
ionizing SED seen by the BLR will then be harder than 
the observed SED, which is likely seen at a smaller inclina- 
tion, closer to a face on view. Assuming the torus-like con- 
figuration where the BLR resides, is restricted to fj, < 0.5, 
and the clear line of sight cone where the AGN is unob- 
scured is restricted to fj, > 0.5, it is plausible to assume 
that the mean inclination of the BLR is fi = 0.3, while our 
line of sight has a mean value of n — 0.8. Table 1 provides 
i'max(/i ~ 0.8)/i^max(/i ~ 0.3), aiid shows that the BLR can 
be exposed to an SED a factor of 1.2 — 2.2 harder than the 
observed SED. If we denote a;o.oi = /iJ'o.oi/fcT'efi.max, and use 
the parameter = a;o.oi(A* = 0.3)/a;o.oi(/i = 0.8), then the 
minimal value for «3ooo derived in eqs.13,14 scales as 
Table 1 lists Af^ values for a range of a, values. Applying this 
inclination correction for the hardness of the ionizing SED 
implies minimal R/3 FWHM values of 9090 km s"^ (a* = 0), 
and 25,820 km s"^ (a* = 0.998). 

The second effect which can make the AD SED harder 
is deviations from the local BB approximation. A deviation 
from a local BB necessarily leads to a photospheric temperar- 
ture T > Teff . Table 2 provides the ratios of I'max /Tiff, max de- 



rived from detailed AD atmosphere models computed with 
TLUSTY. These model are physically equivalent to those 
presented in Hubeny et al. (2000), but utilize the interpola- 
tion scheme described in Davis & Hubeny (2006). Since we 
have a limited range of converged anrmli, we only consider 
TLUSTY models with a. 0.9. In the local BB AD mod- 
els t'max/Tefi.max depends only on a* and //, as M and M 
only shift the SED, and do not affect its shape. In the atmo- 
spheric model the local emission is set by the photospheric 
density and temperature, and Tabic 2 provides models for 
the AD parameters which produce a relatively cold AD, rel- 
evant for our study. Since T > T^s the SED is expected to 
be hotter compared to the BB model for the same parame- 
ters, and thus produce a higher xo.oi value. A notable caveat 
is that some models have a strong Lyman absorption edge, 
which can lead to a lower /ion compared to the local BB AD 
model, despite the higher T. 

Figure 1, lower panel, presents the AD SED derived 
from the TLUSTY models. The deviation of these detailed 
calculations from the simple BB approximation arc strongly 
dependent on the local disk Tgff^niax and the photospheric 
density. The photospheric density also depends (albeit more 
weakly) on the the stress prescription, which is parametrized 
here using a = 0.01 (SS73). The SS73 model gives a verti- 
cally averaged density which scales as a^^ in the radiation 
dominated regime. The dependence of the photospheric den- 
sity on a in TLUSTY models is typically much weaker, for 
reasons discussed in the Appendix of Davis et al. (2006), 
unless a and M are both large enough that the disk be- 
gins to become effectively optically thin. The SED can be 
significantly harder than the local BB model when Tefj^max 
is high. For example, i^max increases by a factor ~ 5 for 
the riiH = 1, Afi — 1 model. However, the significantly 
colder AD which are of interest here, are not far from the 
local BB approximation. Table 2 provides the correction fac- 
tors Aniodci = (TLUSTY) /x(,.oi(BB) for /i = 0.8. Since 
fsooo oc A^''^^ , atmospheric effects can increase the minimal 
R0 FWHM values by a factor of ~ 1.2 - 2. 



3 OBSERVATIONS 

The SED of quasars with a cold AD should show a sharp 
cutoff beyond the Lyman limit. The WLQ population 
comes closest to the lineless quasars population. The WLQ 
with the highest rest frame UV energy probed is SDSS 
J094533.99-I-100950.1, at 2 = 1.66, discovered by Hryniewicz 
et al. (2010). The GALEX FUV photometry extends down 
to rest frame 570A. Figure 2 presents a local BB AD match 
to the overall SED presented by Hryniewicz et al. The 
AD model jjarametcrs were: mg ~ 27, as estimated by 
Hryniewicz et al. based on the Mg II FWHM and the 3000A 
continuum luminosity. We assume an AD with a face on view 
(/i = 1). The value of M was then adapted to reproduce the 
overall SED normalization. The remaining free parameter is 
a,, which is varied to match the FUV turnover. A suitable 
match to the data is found with Mi = 11.4 and o* = 0.3. A 
good match can also be obtained with a higher inclination 
AD model, where the resulting harder SED is compensated 
by allowing a lower a* . A reasonable match can be obtained 
down to fj, = 0.6, and a, = —1, i.e. a BH counter-rotating 
with respect to the AD. A higher inclination, fj, < 0.6, is ex- 
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eluded, as the AD SED becomes too hard, and o* eannot be 
further reduced to compensate for that. Similarly, a value of 
a* > 0.3 is excluded, as the AD is too hard and and ii cannot 
be further increased to compensate for that. These results 
are similar to those obtained by Czcrny ct al. (2011), who 
also discuss the effects of intrinsic reddening on the derived 
AD parameters. We also attempted TLUSTY model fits to 
the data, however the SED of these models is broader than 
the local BB AD models, and a fit of comparable quality to 
the local BB AD models fit could not be obtained. 

Could the observed UV emission be significantly af- 
fected by the foreground Lyman forest absorption? Barger 
& Cowie (2010, Fig. 6 there) show that the cumulative effect 
of the Lyman a absorption systems along the line of sight, 
suppresses the continuum at rest frame 912A< A <1216A by 
a level of 4% for « = 1.1. Using the evolution of the Lyman 
forest systems, dn/dz oc (1 + 2)2-23±i.2 0.9 <« < 1.7 
(Janknecht et al. 2002), the expected continuum suppres- 
sion at z = 1.66 is < 11%, and is thus not significant. 
However, Lyman limit systems can produce significant ab- 
sorption at rest frame A <912A. To affect the observed 
SED the Lyman limit absorber should occur longward of 
observed 1516A (the GALEX FUV effective wavelength), 
i.e have z > 0.66. Using the number density per unit red- 
shift evolution of dn/dz = 0.15 x (1 -|- zY'^ (Songaila & 
Cowie 2010), the expected number of absorbers integrated 
over 0.66 < z < 1.66, is 0.66. A Lyman limit absorber will 
likely cause a break in the spectral curvature, which is not 
seen in the data (Fig. 2). A further complication may arise 
if a foreground Lyman break is diluted by a Lyman contin- 
uum emission edge produced by the AD. UV spectroscopy 
is required to clearly exclude foreground absorption effects 
on the observed FUV SED of SDSS J094533.99-M00950.1. 



4 DISCUSSION 

The clear prediction of local BB AD models is that AGN 
with the correct combination of M and I harbor a lumi- 
nous, but non-ionizing, continuum source. The UV SED of 
an AGN with /ion = 0.01 will show a peak at A > 3200A, 
and a steep drop at A < lOOOA, and may form a lineless 
quasar, depending on the relative strength of the X-ray ion- 
izing power-law continuum. An AGN with /ion = 0.1 will 
show a peak at A > I8OOA, and will likely form a WLQ, 
with an SED similar to that observed in the WLQ SDSS 
J094533.99-I-100950.1. 

A cold non-ionizing AD is expected in luminous, Lboi ^ 
lO''^ erg s~^, quasars even at relatively high values of I < 0.2, 
in particular for a, ~ 0. In low luminosity AGN, Lboi 
10*"^ erg s~^, the SED is predicted to become non- ionizing 
only at very low values of / < 10~^. Since the X-ray/optical 
luminosity ratio is observed to increase with decreasing Lboi, 
such low luminosity and low / AGN may be characterized 
by more prominent line emission, excited by the relatively 
stronger X-ray power-law emission. 

What is the nature of the expected BLR emission? The 
flat X-ray power-law continuum produces an extended re- 
gion of low ionization in the BLR gas, which is expected 
to cool mostly through coUisionally excited low ionization 
lines, such as Fe II, Mg II, and Balmer lines (coUisionally 
excited due to highly trapped Lya photons). The narrow 



line region gas is expected to cool mostly through forbidden 
low ionization lines, such as O I, N II, and S II lines, as seen 
in LINERs (see photoionization model results in Ferland & 
Netzer 1983). 

4.1 Are we missing lineless quasars? 

Is it possible that AGN surveys are missing a significant 
number of moderately high high M AGN which are line- 
less?. This can be addressed by surveys of AGN based on 
color, independent of line emission. This was carried out 
by Diamond-Staruc (2009) using the SDSS sample, where 
they defined WLQ as having a Lya EW<15.4A, which is 3(7 
from the mean Lya EW=63.6A. They find that the fraction 
of quasars classified as WLQ increases from 1.3% at « < 4.2 
to 6.2% at > 4.2. Thus, lineless quasars do not form a 
significant part of the quasar population, as expected given 
the required very high M values in luminous AGN, in order 
to have a cold AD. 

4.2 Do LINERs have a cold AD? 

In contrast with luminous AGN, a significant fraction of low 
luminosity AGN near the more massive black hole should fall 
in the cold AD regime. Theoretically, it is not clear whether 
the thin AD configuration remains valid at Z < 10~^, or 
whether the accretion transforms to an optically thin and 
geometrically thick configuration, which may produce a hard 
power- law continuum source (e.g. Narayan et al. 1998). Ob- 
servationally, the narrow emission lines of low luminosity 
AGN in massive galaxies are generally characterized as LIN- 
ERs, and are likely to be powered by a hard continuum 
source (Ferland & Netzer 1983), as their observed SED 
apparently suggests (e.g. Ho 2008, Eracleous et al. 2010a; 
cf. Maoz 2007). Thus, the significant contribution of X-ray 
heating in low luminosity AGN, given their flat aox, may 
imply that FUV heating is not significant in these objects 
anyhow, and its absence as the AD becomes cold will not af- 
fect their line emission significantly, unlike typical luminous 
AGN where the FUV dominates the line excitation. The 
line emission in LINERs may also be powered by additional 
mechanisms (e.g. Eracleous et al. 2010b), which further low- 
ers the diagnostic power of the line emission on the ionizing 
SED shape. 

Interestingly, Lawrence (2005) finds that the change in 
the SED of a few potentially characteristic AGN, extend- 
ing over a range of 10® in luminosity, is consistent with the 
expected change for an optically thick AD at a fixed M. 
It remains to be explored whether the optical-UV SED in 
the lowest luminosity AGN (e.g. LINER) can be fit by a 
cold AD, or whether only a featureless hard power-law is re- 
quired. A cold AD in low luminosity AGN cannot produce 
a lineless AGN, unless Oox is unusually steep, and ionizing 
radiation is indeed missing. 

4.3 True type 2 AGN 

In "true type 2 AGN" which show normal (non LINER) 
narrow lines, the absence of broad lines does not result from 
a cold AD, as an apparently normal photoionizing radiation, 
which excites the narrow lines, is present. In such objects. 
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the BLR may be absent due to the absence of the mechanism 
which produces it, such as an AD wind, as proposed by 
Nicastro (2000), Nicastro et al. (2003), Czerny et al.(2004), 
EUtzur & Shlosman (2006) and EUtzur & Ho (2009). 

4.4 What produces WLQ? 

Are most WLQ powered by a cold AD? The archetype 
WLQ, PG 1407+265, is probably not, as it's optical-UV 
SED appears to bo similar to the average quasar SED (Mc- 
Dowell et al. 1995). It does not show the expected peak at 
A > I8OOA and a sharp drop at A < lOOOA (although the 
the S/N at A < lOOOA is rather low). However, its peculiar 
emission line properties, with some detectable low ioniza- 
tion lines (Balmer lines, Mg II, Fe II), but very weak or no 
detectable higher ionization lines, is characteristic of other 
WLQ (Leighly et al. 2007a; Diamond-Stanic et al. 2009; 
Plotkin 2010b; Shcmmcr ct al. 2010; Hrynicwicz et al. 2010; 
Wu et al. 2011). In the case of PHL 1811, this line emission 
is modeled by an unusually soft ionizing spectrum (Leighly 
et al. 2007b), which is however qualitatively similar to the 
emission of a cold AD. The soft ionizing SED produces little 
heating per ionization, and thus weak coUisionally excited 
lines compared to recombination lines. 

Some WLQ may be driven by a different processes. The 
C IV EW is tightly inversely related with I (Baskin & Laor 
2004), and also with other emission line properties, such as 
the [O 111] and Fc II EW (Baskin & Laor 2005), which 
are part of the eigenvector 1 correlations of various emission 
line properties (Boroson & Green 1992). The G IV EW is 
also related to the C IV line peak velocity shift and profile 
asyrmnotry (e.g. Baskin & Laor 2005; Richards ct al. 2011), 
and aox (Baskin & Laor 2005; Wu et al. 2009). This large 
set of correlations suggests that local effects in the BLR (e.g. 
metallicity, kinematics), may control the EW of C IV, and 
possibly of other high ionization UV lines. If the driving 
mechanism for WLQ is a filtered ionizing continuum which 
illuminates the BLR (e.g. Richards et al. 2011; Wu et al. 
2011), then the effect on the BLR excitation may be similar 
to the effect of a cold AD. However, the trend of decreasing 
C IV EW with increasing / cannot be confused with the cold 
AD effect, as higher I implies a hotter AD SED, rather than 
a colder one. 

4.5 The absence of AGN with very broad lines 

The combination of M and M which produces a cold AD, 
can be transformed to a limit on the H/3 FWHM, with 
very weak Lopt dependence. The local BB AD may be- 
come cold at V > 8,000 km s~^ when a, = 0, and for 

V > 16,000 km s"^ when a, = 0.998. The lower limit on 

V can be higher by a factor of up to 2 due to the inclination 
dependence of the ionizing continuum, and possible devia- 
tion of the disk emission from a local BB continuum. AGN 
do show a sharp drop in their rmmber at t; > 10, 000 km s~^ 
(e.g. Hao et al. 2005, Shcn et al. 2008), and this drop is 
independent of luminosity (Shen et al. 2008, Fig. 3 there). 
Laor (2003) suggested there may be a yet unknown physical 
mechanism which suppresses the BLR at > 25, 000 km s^^ . 
The AD becoming too cold to ionize the BLR, may provide 
a natural explanation for the observed luminosity indepen- 
dent limiting v, beyond which broad lines are not seen. 



How does the line emission of AGN with the broad- 
est lines appear? Do such objects enter the WLQ regime? 
Such a systematic study is not available yet. However, very 
broad line objects tend to have a double-peaked emission 
profile, and the spectrum of the "prototypical" such object, 
ARP 102B, shows strong Balmer and Mg II lines, strong low 
ionization forbidden lines, and relatively weak higher ion- 
ization lines (Halpern et al. 1996). Only anecdotal informar 
tion is available for a couple of more such objects (Storchi- 
Bcrgmann et al. 2005; Eracleous ct al. 2009), partly con- 
sistent with the results for ARP 102B. A systematic study 
is clearly required to explore whether very broad line AGN 
indeed have emission line properties intermediate between 
normal AGN and WLQ. 

What is the SED of very broad line quasars? Strateva 
et al (2008) shows the overall SED of five double-peaked 
Balmer lines AGN, with v ~ 14, 000 - 20, 000 km s"^ Their 
optical-UV SED is generally consistent with the average 
AGN SED, but the FUV at A < lOOOA is not probed. The 
objects have a flatter aox than the average for quasars at a 
similar luminosity (see also Strateva et al. 2006). This may 
indicate indirectly that significant X-ray heating is required 
for generating the broad lines in these objects, as very broad 
line quasars with an average aox may be selected against due 
to their weaker line emission. The flatter Oox may therefore 
compensate for a smaller contribution of the UV continuum, 
if their AD is indeed colder. A systematic study of the FUV 
emission of very broad line AGN can test if their FUV is 
indeed steeper. 



4.6 Continuum based searches for lineless quasars 

An alternative approach to test the cold AD scenario is to 
explore the emission line properties of AGN where the SED 
indicates weak FUV emission. It is important to note that 
red quasars (e.g. Richards et al. 2003) are not suitable ob- 
jects, as a cold AD produces a normal blue continuum at 
A > 32OOA. Red quasars are either affected by dust ex- 
tinction, or by some other modification of the continuum 
emission mechanism. Gold AD candidates can be found by 
looking for quasaxs with a blue NUV slope, a ° > — 1, 

A ^3200 A 

and a red FUV slope, a ° < —2, as seen in SDSS 

^ A<1000A 

J094533.99+100950.1 (Fig.2). Such observation cannot be 
achieved from the ground, as the required z > 3 to probe the 
FUV slope from the ground, leads to a high probability for 
an intervening Lyman limit absorption system (dn/dz > 2). 
Space based UV observation of ^ < 2 quasars are essen- 
tial for that, as was done by Hryniewicz et al. (2010) using 
GALEX. 

The observed optical-UV SED of optically selected 
quasars shows a small dispersion (Elvis et al. 1994; Richards 
et al. 2006). Due to this Davis & Laor (2011) deduced that 
the accretion efficiency, rj, increases with M, to compensate 
for the drop in the FUV with increasing M for AD mod- 
els at a fixed r/. One may worry that the implied rise in 
r) with M is just a selection effect, as high M and low 77 
quasars necessarily have a cold AD, produce weak line emis- 
sion, and arc selected against in a broad emission line AGN 
survey. However, a color based AGN survey is not blind to a 
cold AD quasar, and as noted above, Diamond-Stanic (2009) 
find that only 1.3% oi z < 4.2 color selected AGN are WLQ. 
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Thus, there is no large population of low a« high M quasars, 
and the high mean rj (high a*) values for high AI AGN de- 
duced by Davis & Laor (2011), is not driven by emission line 
selection. 

Wc note that quasars generally show a soft-excess X-ray 
emission below 1 keV, compared to the 2-10 keV power-law 
emission (Wilkes & Elvis 1987), which may extend as a sin- 
gle power-law component to the FUV (Laor et al. 1997). 
This component is commonly attributed to comptonization 
in a warm surface layer above the AD (e.g. Kriss ct al. 1999). 
Such a layer can turn a cold AD SED into an ionizing SED. 
The SED of SDSS J094533.99-I-100950.1 is consistent with 
the exponential cold AD FUV cutoff, and shows no indica- 
tion for a power-law component. Such a warm comptonizing 
surface corona may be lacking in a cold AD (Czerny et al. 
2011). 

4.7 Cold AD in Blazars 

Ghisellini et al. (2009, 2010) attribute the optical-UV emis- 
sion feature in blazars observed with the Swift satellite to 
thermal AD emission, based on the low polarization and 
low variability of this component. This component shows 
a peak in the NUV and a steep falloff in the FUV, and is 
fit with an AD model with very high values of M, up to 
in some objects. These are therefore generally 
cold AD models, and if the BLR emission in these objects 
is powered mostly by the unbeamed thermal disk emission, 
some of these blazars should appear as WLQ. Indeed, some 
of the objects in Ghisellini et al. (2010) do show weaker UV 
hnes (e.g. S5 0014-1-81, see Kuhr et al. 1983, Sargent et al. 
1989; PKS 2149-306, see Wilkes 1986; RBS 315, see Elli- 
son et al. 2008), which may provide independent support 
for the cold AD interpretation of their opticaJ-UV bump. 
However, intervening Lyman limit absorption systems are 
clearly contributing to the observed FUV steep falloff in 
these z ~ 2.5 — 3.5 blazars, so the intrinsic AD continuum 
may not be as cold as observed. In addition, one cannot 
exclude additional line excitation from the highly luminous 
jets in these objects. 



powered by a cold AD. The best strategy is to obtain UV 
spectra of z ~ 1 — 2 very broad line quasars, to minimize 
the likelihood of intervening Lyman limit absorption by the 
IGM, or to detect such absorption if it occurs. 

Color selected AGN samples indicate that only a small 
fraction (< 6%) of luminous AGN may harbor a cold AD. 
This can result from the high mean ij values for AGN at high 
M, which may imply a higher o* and a hotter AD emission 
which produces a higher /ion- 

Searches for quasars with a blue optical-NUV and a red 
FUV continuum, using ground based optical spectra, and 
space based UV spectra of 1 < 2 < 2 quasars, may be an 
efScicnt mechanism to reveal cold AD AGN. The continuum 
of a WLQ provided by Hryniewicz ct al. (2010), which ex- 
tends down to A = 570A, is surprisingly well fit by a simple 
local BE AD models. This is in contrast with most AGN, 
where the SED cannot be well fit by a simple local BB AD 
model. The local BB models provide a better match than 
the TLUSTY-based models (Hubcny et al. 2000), suggest- 
ing that the effects of the Lyman edge and electron scatter- 
ing opacity on the spectrum may be overestimated by these 
models. 

Since both M and Lopt are determined from the obser- 
vations of SDSS J094533.994-100950.1, only a, and ^^ remain 
as free parameters to match the observed SED. The fit allows 
us to place an upper limit on the black hole spin (a* ^ 0.3), 
and a lower limit on the inclination of the AD (n > 0.6). 
The precise limits are somewhat uncertain, in part because 
the Mg II line-based estimates of M are uncertain. A larger 
M would allow a wider range of o, and fi, but a smaller M 
would provide tighter constraints: am even lower limit on a* 
and higher limit on /x, as further discussed by Czerny et al. 
(2011). 

It is worthwhile to explore whether there are additional 
objects which can be well fit by a cold AD SED, with only 

a, and fi as free parameters. If such objects exist, the SED 
may provide a useful tool to constrain a, and fi in such AGN, 
as commonly done in X-ray binaries. The AD fit may also 
provide some new insight on the structure of the accretion 
disk atmosphere. 



5 CONCLUSIONS 

Standard thin AD models imply a non ionizing continuum 
for high M (> 3 X IO'^Mq), even for fairly luminous quasars 
(iopt = 10"*^ erg s~^), in particular if the black hole is 
non-rotating. Such cold AD may reside in AGN where the 
Balmer line FWHM> 8,000 - 24,000 km s"^ The ioniz- 
ing continuum in such objects originates only in the X-ray 
power-law component, which will likely produce an extended 
partially ionized region in the illuminated gas, which cools 
mostly through low ionization lines. If < O.lLboi, as com- 
monly seen in luminous AGN, the ionizing luminosity will 
be weaker by a factor of > 5 than in average AGN. Such 
objects will appear as WLQ. If « Lboi such quasars 
may be lineless. 

Cold AD provide a natural explanation for the sharp 
drop in the number of AGN with broad line FWHM> 
10, 000 km s~^, which is independent of luminosity. FoUowup 
FUV observation axe required to test if very broad line AGN 
show the expected spectral turnover in the FUV, if they are 
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Figure 1. The SED of various AD models. All models have M\ = 1 and ^ = 0.8. The vertical dashed line shovfs the threshold v for 
ionizing H. The two upper panels show various local BB accretion disk models for a» = and a* = 0.998. The SED gets softer as M 
increases, and as a, decreases. The lower panel shows the TLUSTY model (Hubeny et al. 2000). with a = 0.01. The SED is harder than 
the local BB model, and the effect is larger when Tgfj is larger. The fraction of ionizing radiation, /iom is raised by the atmospheric 
effects, compared to the local BB AD models. 
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Figure 2. A local BB AD model match to the z = 1.66 weak emission line quasar SDSS J094533. 99+100950.1 from Hryniewicz et al. 
(2010). The TLUSTY models produced a lower quality match to the data. The photometric data points are taken from non simultaneous 
observations by 2MASS (triangles), SDSS (squares), and GALEX (diamonds). The AD model fit parameters are listed. The value of M 
is taken from Hryniewicz et al. The model peak position is mostly set by the combination of // and a* , and the normalization mostly by 
M. The AD model fit gives /ion = 0.06, above our definition of /ion = 0.01 for potentially lineless quasar, and this higher value may be 
driving the weak low ionization line emission observed. 



Table 1. 



a* 


^eff,max/2o 


^^'max/A^^eff ,niax 


fO.l/f'max 


I'O.Ol/t'max 


t'max(Al = 0.8)/fc'max(Al = 0.3) 







0.1 


2.18 


1.91 


3.43 


0.85 


1.17 


0.9 


0.23 


1.77 


1.97 


3.54 


0.64 


1.59 


0.998 


0.45 


1.18 


2.03 


3.54 


0.46 


2.03 



The results for local BB AD models, observed at ;U = 0.8. All the results in the table are independent of the 
values of mg and Mi. 
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Table 2. 



a* ms Ml /ll'O.l/feTefl^niax hl/O.Ol/kTsfl^ialiyc ^model 






10 


1.22 


6.1 


9.9 


1.32 





32 


12.2 


6.8 


11.6 


1.56 





100 


123 


7.2 


16.0 


2.15 


0.9 


10 


0.28 


6.6 


12.1 


1.93 


0.9 


32 


2.84 


8.3 


14.5 


2.32 


0.9 


100 


28.4 


8.9 


16.6 


2.67 



The TLUSTY model results, observed at = 0.8. All models have 
a fixed nifi X Mi = 1, and therefore a fixed Tq. For the BB model, 
h;/0.i/fcTcfi,max = 4.2 - 3.5, and /i!^o.Ol/fcTeff,max = 7.5 - 6.3, for 
a» = 0—0.9. In contrast to the BB case, we do not compute a, = 0.998 
models because our TLUSTY model atmosphere table is not extensive 
enough to reliably compute disk spectra for such high spins. 



